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Abstract 

We analyze the quantum-mechanical rotational excitation/de-excitation spectrum and cross sec- 
tions of CN molecules during low and high-energy collisions with protons, H + . The problem is of 
significant importance in astrophysics of the early Universe, specifically connected with the prob- 
lems of cosmic microwave background (CMB) radiation. A quantum-mechanical close-coupling 
Q-i method is applied in this work. The cyanide molecule (CN) is treated as a rigid rotor, i.e. the 

distance between the carbon and nitrogen atoms is fixed at an average equilibrium value. The new 
results of the excitation/de-excitation cross-sections and corresponding thermal rate coefficients are 
compared with the results of few previous calculations performed on the basis of few approximate 
semiclassical frameworks. The interaction potential between CN and H + is taken in the follow- 
ing form: proton induced polarization potential + proton-dipole potential + proton-quadrupole 
potential. 
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I. INTRODUCTION 



Matter emits electromagnetic radiation when it has a temperature above zero Kelvin. 
This is also called thermal radiation because it is generated by thermal motion of particles 
in the matter. Electromagnetic radiation can be absorbed by matter as well and matter 
that absorbs all radiation is called a blackbody. A blackbody has a uniform temperature 
and produces radiation characterized by the frequency distribution that depends on its 
temperature. Though a background radiation, referred to as excess noise, was noticed as 
early as 1941, it wasn't until 1965 that Arno Penzias and Robert Wilson (5)] used highly 
sensitive equipment and careful observation, which led to the labeling of this radiation as 
cosmic microwave background (CMB) radiation: a nearly isotropic (even in all directions) 
thermal radiation in the Universe. It is know that when a molecule absorbs a photon it 
moves from a ground state to an excited energy state. Molecules in the space between 
Earth and distant stars produce absorption lines seen on Earth that correspond to the 
spectra of those stars. Most of these lines are viewed at the ground state of the molecules, 
however, researchers found that CN molecules observed in various parts in space all have faint 
absorption lines at the first excited energy state 4J . Questions arose as to what could excite 
the CN molecules in space? The temperature of radiation can be found from the amount of 
molecules excited by that radiation. Measurements of the spectrum of this radiation showed 



a. 



that it was almost identical to blackbody radiation with the temperature of about 3K 
Since these measurements where observed for almost all CN molecules in different parts of 
the Universe, it showed that whatever was exciting the molecule was distributed fairly evenly 
throughout space. This radiation cannot be attributed to energy from stars or intergalactic 
gas, for example, since this radiation's spectrum does not look like the spectrum for stars 
or intergalactic gas or etc. And also, if the energy of this radiation across the Universe is to 
be summed up, it would be much too great to be radiation from stars or such. Therefore, 
researchers believe that whatever is exciting the CN molecules in space ,with a spectrum of 
a blackbody at 3K, is CMB. 

As mentioned previously, CMB is isotropic, however, not perfectly so. It has an anisotropy 
of about 0.1 percent, which is due to a Doppler shift caused by the motion of the solar 
system through CMB. This anisotropy is very small, therefore, CMB is still regarded as 
almost perfectly uniform radiation in the Universe. One might wonder what caused CMB 



radiation to occur in the Universe at all? To answer this question, researchers have focused 
on the blackbody characteristic of the radiation, since this is an intrinsic characteristic. 
Knowing the present temperature of the radiation being at about 3K, one can calculate the 
temperature at an earlier time, which turns out to be quite large. This large temperature 
greatly points to the existence of the the Big Bang: the rapid expansion of the once a 
dense and compact Universe. Another model that goes together with the Big Bang theory 
is the model of the expanding Universe. It is strongly shown experimentally and supported 
theoretically that the Universe is expanding. In the 1920's Edward Hubble found that 
spectra of distant galaxies are redshifted, which means the galaxies are traveling away from 
Earth. The model of the expanding Universe states that the cold Universe of today must have 
had very high temperature and much higher density when it first begun billions of years ago. 
With temperatures greater than 10 8 K, the Universe was once a plasma at thermodynamical 
equilibrium and made of electrons, protons, and helium ions, that radiated and absorbed 
photons. As the Universe expanded, the temperature of the plasma and radiation decreased. 
At temperatures of about 4000K, recombination of protons and electrons occurred and the 
plasma became a mixture of neutral hydrogen and helium atoms. The Universe became 
transparent to radiation, as it is now. As the Universe continued to expand, the temperature 
of the radiation declined but the blackbody character of the radiation remained as a "relic." 
CMB radiation is presumed to be "relic" radiation from the early times of the Universe. 

II. QUANTUM-MECHANICAL APPROACH 
A. Equations and Methods 

A quantum mechanical method was used in deriving our results for rotational quantum 
transitions in collisions between the CN molecule and a proton H + , that is: 

CN(j) + H + ^H + + CN(j'). (1) 

The interaction of H + + CN is treated as a rigid rotor model. We used close-coupling 
quantum-mechanical methods to calculate the cross sections and collision rates of a CN 
molecule with a proton H + . The Schrodinger equation for an a + be collision in the center 
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of mass frame, where a (H + ) is an atom and be (CN) is a linear rigid rotor, is: 

where Pg is the relative momentum between a and be, Mr is the reduced mass of the atom- 
molecule (rigid rotor in this model) system a + be: Mr = m a (m^ + m c ) / (m a + mj + m c ), \x 
is the reduced mass of the target: fi = m^nicj (777,6 + mc), ? is the angle of orientation of the 
rotor ab, V(r, R) is the PES for the three-atom system abc, and E is the total energy of the 
system. 

To solve Eq.(2), the following 3-atom wave expansions are used [111 ]: 



*M)= E ^^^M), (3) 

JAfjL 



and 
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r, fl) = ^ CjwLLmz Yjmi (r)Y Lm2 (R) , (4) 

A set of coupled differential equations are achieved by solving the Schrdinger equation 
through the use of the free software MOLSCAT (ver.14, J. Hudson and S. Green (1994)): 

^ - ^jP 1 + k%) U J 3 t\R) = 2M R g / < 0/f(f, R)\V(r, R)\<f>™(r, R) > 

xU^ L I ,(R)dfdR. (5) 

The formulas used for analyzing the cross sections are related to the asymptotic behavior 
of the H + + CN wave function using U/ T M (R — > 00). The numerical results are fitted to the 

n 

known asymptotic behavior of Uj^(R) from [121 ] : 



and 



*-("/', J» = * J] (2 J + - S J (j'L';jL; E)\ 2 . (7) 
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The kinetic energy is e = i? — B e j(j + 1), where t3 6 is the rotation constant of the rigid rotor 
be, i.e. the hydrogen molecule. 



B. CN + H + Interaction Potential 

The PES equations used were found from works Where U p is the polarization 



potential, Ud is the proton-dipole potential, and U q is the proton-quadrupole potentia 
values for the potential parameters for these equations were found and referenced in [2 



. The 



U p+C n(7) = U p {f) + U d {f) + U q {f) (8) 



-e 2 a ■ (r 2 
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/2, if r < r s 



U P (r) = I -e 2 a /2(r 2 + r 2 ) 2 - [e 2 a ■ (r - r s f ■ P 2 (cos(9))]/[2(r 2 + r 
(b 2 + (r - r s ) 2 )], if r > r 8 , 
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UJr) 



0, if r < r s 

(eD ■ (r - r.fP^cos^)))/^ 2 ■ [b 2 + (r - r s ) 2 ]), if r > r s 



0, if r < r s 

eQ/r 3 ■ (r — r q ) 2 P 2 (cos(6))(b 2 + (r — r 9 ) 2 ) -1 , if r > r s , 

Where r is the vector that joins the center of mass of the CN molecule to the H + , and 9 is 
the angle between the r and the axis of the CN molecule. The potential surface is modeled 
in Fig. 1. 



III. RESULT AND DISCUSSION 



We have taken test calculations to insure that the convergence of the results with respect 
to all parameters that enter in the propagation of the Schrodinger equation exists. The 
parameters include the atomic-molecular distance R, the total angular momentum J, the 
number of total rotational levels to be included in the close-coupling expansion and others. 
Attention has been given to the total number of numerical steps in the propagation over 
the distance R of the Schrodinger Eq. flS]). The parameter R ranges from 0.75 a.u. to 
20.0 a.u. Up to 50000 propagation points have been used as well as different mathematical 
propagation schemes included in MOLSCAT. 
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Fig. 2 shows our results for the rotational energy transfer cross-sections of the process in 
Eq. (Tjjj when j transitions from to 5. The results converge better at higher total energy. 
Fig. 3 shows separately the results for j transitions from to 1 as compared to work [1]. We 
can see from Fig. 3 that the pattern of the graph obtained from our results is similar to 
work [1], but not exact. In Fig. 4 and Fig. 5, we compare different J values from our data. 
From the graphs, we can see that the cross sections are close. 

In summary, the subject of this work is the quantum-mechanical study of the state- 
resolved rotational relaxation and excitation cross sections in ultra-cold collisions between 
CN molecules and protons H + . We analyzed the excitation of the CN molecule. A model 
3-body PES for H + -CN has been adopted from a previous work [l| and the potential is 
shown in Fig. 1. Calculation for total elastic scattering cross sections and for low quantum 
rotational transition states have been performed. A test of the numerical convergence was 
taken. 

In the future, we plan to compute the thermal rate coefficients. 
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FIG. 1. The H + + CN potential energy surface used in this work. 
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FIG. 2. The H + + CN Rotational Energy Transfer Cross-sections 
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FIG. 3. Results for the H+ + CN (0) to H+ + CN(1) 
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<D— 1> Result with J=2000 
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FIG. 4. Results for the H+ + CN (0) to H+ + CN(2) 
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FIG. 5. Results for the H+ + CN (0) to H+ + CN(4) 
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